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PROCEDE ET DISPOSITIF D'IMAGERIE UTILISANT LES ONDES DE CISAILLEMENT. 

Proc6d6 d'imagerie pour observer la propagation 
dljne onde impulsionnelle de cisaillement basse frequence 
simultan£ment en une multitude de points d'un milieu vis- 
coeiastique diffusant (1). A cet effet, on 6met k cadence ul- 
trarapide des ondes ultrasonores de compression qui 
permettent d'obtenir une succession d'images du milieu, 
puis on traite en temps differ^ les images ainsi obtenues par 
intercorrGlation, pour determiner en chaque point de chaque 
image les mouvements du milieu lors de la propagation de 
I'onde de cisaillement. 
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Prncede et disposifcif d'imaqerie utilisant — — <?ttd<?S fle 
r r iaaillement. 

La presente invention est relative aux procedes et 
5 dispositifs d'imagerie utilisant les ondes de cisaille- 
ment . 

Plus particulidrement , 1' invention concerne un 
procede d'imagerie utilisant les ondes de cisaillement 
pour observer un milieu viscoelast ique diffusant qui con- 
10 tient des particules ref lfichissant Tes" ondes ultrasonores 
de compression, procede dans lequel on genere une onde 
elastique de cisaillement dans le milieu viscoelast ique et 
on observe, au moyen d 1 au moins une onde ukrasonore de 
compression, le deplacement du milieu viscoelast ique sou- 
15 mis a ladite onde de cisaillement. 

Le document US-A-5 810 731 decrit un exemple d'un 
tel procede, dans lequel 1 ' onde de cisaillement est gene- 
ree localement a 1 • interieur du milieu viscoelast ique ob- 
serve, au moyen de la pression de radiation d'une onde ul- 
20 trasonore modulee et focalisee sur un point a observer. On 
envoie ensuite sur ce point focal une onde ultrasonore 
supplementaire dont la reflexion permet de connaitre cer- 
tains parametres de propagation de 1 • onde de cisaillement 
(notamment la viscosite dynamique du milieu et son module 
25 de cisaillement) au niveau du point focal susment ionne . 

Cette technique presente 1 1 inconvenient de permet - 
tre 1 ' analyse d'un seul point du milieu viscoelast ique 
etudie a chaque fois qu'une onde de cisaillement est gene- 
ree. Si 1 1 on veut obtenir une image complete du milieu 
30 viscoelastique observe, il est necessaire de repeter 
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l'operation un tres grand nombre de fois, ce qui implique 
un temps de pause important (par exemple, plusieurs minu- 
tes) pour obtenir cette image. 

Ce temps de pause important rend ce procede de 
1 'art anterieur peu pratique a utiliser. 

De plus, un tel temps de pause peut nuire a l'uti- 
lisacion dudit proced6 pour obtenir une image d 1 un tissu 
vivant, qui est tou jours en mouvement . 

La presente invention a notamment pour* but 'de pal- 
lier ces inconvenients. 

A cet effet, selon 1' invention, un procede du 
genre en question est essent iellement caracterise en ce 
qu'on genere 1 1 onde de cisaillement en ; appliquant au "mi- 
lieu viscoelastique une excitation ayant ' la forme d'une 
impulsion basse frequence qui presente une frequence cen- 
trale f comprise entre 20 et 5000 Hz, cette impulsion 
basse frequence pr£sentant une duree comprise entre 1/2 f 
et 20/f, 

en ce qu'il comporte une etape d 1 observat ion de propaga- 
tion au cours de laquelle on observe la propagation de 
1 1 onde de cisaillement simultanement en une multitude de 
points dans le milieu observe, ces points formant un champ 
d ' observat ion sensiblement continu s 1 etendant au moins se- 
lon un premier axe, cette etape d 1 observation de propaga- 
tion de 1 ' onde de cisaillement consistant a : 

emettre dans le milieu observe une succession 
d'au moins 10 tirs d'ondes ultrasonores de compression a 
une cadence comprise entre 100 et 100 000 tirs par se- 
conde , 
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detecter et enregistrer en temps reel les echos 
generes par les particules ref lechissantes du milieu vis- 
coelastique a chaque tir d'onde uitrasonore, ces echos 
correspondant (directement ou indirectement) a des images 
successives du milieu observe, 

et en ce que ledit procede comporte en outre une etape ul- 
terieure de traitement d' image au cours de laquelle on 
craice en temps differe les images ainsi obtenues au moins 
par intercorrelat ion entre images successives, pour deter- 
miner eh "cfiaque point du champ d' observation un parametre 
de mouvement choisi entre le deplacement et la deformation 
du milieu viscoelastique , de fagon a obtenir ainsi une 
succession d' images montrant 1' evolution du parametre de 
mouvement du milieu viscoelastique sous l'effet de la pro- 
pagation de l'onde de cisaillement . 

Grace a ces dispositions, on obtient un film il- 
lustrant clairement la propagation de l'onde de cisaille- 
ment dans le milieu viscoelastique, qui peut permettre par 
exemple, dans des applications medicales, de reperer di- 
rectement des zones cancereuses dans les tissus d'un pa- 
tient : la propagation des ondes de cisaillement s'y de- 
roule en effet tres differemment des zones voisines . 

Ce reperage s'effectue beaucoup plus facilement 
que par une observation classique par simple echographie 
uitrasonore, puisque la propagation des ondes de cisaille- 
ment est fonction du module de cisaillement du milieu, 
lui-meme tres variable entre une zone de tissus sains et 
une zone de tissus cancereux : le module de cisaillement 
varie typiquement dans un rapport de 1 a 30 entre une zone 
saine et une zone cancereuse, alors que le module de com- 
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pression, qui regit la propagation des ondes acoustiques 
de compression utilisees dans 1 ■ echographie ultrasonore, 
varie seulement de 1 1 ordre de 5% entre un tissu sain et un 
tissu cancereux. 

On notera que le film obtenu illustre la propaga- 
tion de 1 ' onde de cisaillement beaucoup plus clairement 
que la simple succession des images donnees par les parti- 
cules ref lechissantes du milieu, puisque ledit film permet 
de visualiser a chaque instant les zones du milieu observe 
qui subissent des mouvements de mime ampleur du fait de la 
propagation de 1 1 onde de cisaillement/ alors que la suc- 
cession des images des particules ref lechissantes ne per- 
mettrait de visualiser qu'un brouillard de points lumineux 
en mouvement . 

Dans des modes de realisation preferes du procede 
selon 1' invention, on peut eventuellement avoir recours en 
outre a 1 1 une et/ou a 1' autre des dispositions suivantes : 

la duree de 1' impulsion basse frequence est 
comprise entre l/2f et 2f ; 

la frequence centrale de 1' impulsion basse fre- 
quence est comprise entre 30 et 1000 Hz / 

le milieu viscoelast ique observe est constitue 
par un corps vivant comprenant au moins un organe interne 
soumis a des mouvements impulsionnels , 1' impulsion basse 
frequence qui genlre 1 ' onde de cisaillement etant consti- 
tute par un mouvement impulsionnel dudit organe interne ; 

le milieu viscoelast ique observe est delimite 
par une surface exterieure, et 1' impulsion basse frequence 
est appliquee au niveau de cette surface exterieure ; 
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1' impulsion basse frequence est appliquee par 
un moyen d* excitation choisi parmi : 

une onde acoustique generee par au moins un trans - 
ducteur acoustique, 
5 et un choc gen§re localement par contact physique 

au niveau de la surface exterieure du milieu viscoelasti- 
que ; 

on emet les tirs d'ondes ultrasonores de com- 
pression et on detecte les echos generes par les particu- 

10 les ref lechissantes du milieu viscoelast ique , au moyen 
d'une batterie de transducteurs qui comprend au moins un 
transducteur et qui est disposee au contact de la surface 
exterieure du milieu viscoelast ique , 1 ' onde de cisaille- 
ment etant appliquee au milieu viscoelast ique en imposant 

15 un deplacement impulsionnel a ladite batterie de transduc- 
teurs ; 

ledit parametre de mouvement est la deformation 
du milieu viscoelastique : cette disposition est particu- 
lierement utile dans le dernier cas envisage ci-dessus, 
20 puisqu'elle permet de s'affranchir du deplacement de la 
batterie de transducteurs, deplacement qui perturberait 
sinon la mesure du deplacement des points du champ d' ob- 
servation ; 

au cours de 1 1 observat ion de la propagation de 
25 l'onde de cisaillement , on emet entre 100 et 10 000 tirs 
d'ondes ultrasonores de compression a une cadence comprise 
entre 100 et 100 000 tirs par seconde ; 

le champ d 1 observation s'etend au moins selon 
un plan comprenant d'une part, le premier axe et d f autre 
30 part, un deuxieme axe perpendiculaire au premier axe ; 
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au cours de l'etape d ' observation de propaga- 
tion, on utilise une batterie de piusieurs transduc teurs 
acoustiques disposes au moins selon le deuxieme axe pour 
emetrre les tirs d'ondes ultrasonores de compression et 
deteccer les echos generes par les particules ref lechis - 
santes du milieu viscoelast ique , les echos detectes par 
chaque transducteur acoustique etant memorises directement 
sans traitement prealable au cours de l'etape d 1 observa- 
tion de propagation, et l'etape de traitement d'image com- 
prenant une sous -etape preliminaire de formation de voies 
au cours de laquelle on genere une image du milieu visco- 
elastique correspondant a chaque tir d ' onde ultrasonore.de 
compression, par combinaison d'au moins certains des echos 
re^us par les differents transducteurs ; 

l'etape de traitement d' images est suivie (im- 
mediatement ou non) par une etape de visualisation au 
cours de laquelle on visualise au ralenti un film consti- 
tue par la succession des images traitees, chaque point de 
chaque image presentant un parametre opt ique qui varie 
suivant la valeur du parametre de mouvement affecte a ce 
point ; 

le parametre optique est choisi parmi le niveau 
de gris et le niveau chromatique ; 

l'etape de traitement d' images est suivie (im- 
mediatement ou non) par une etape de cartographie au cours 
de laquelle, a partir de 1' evolution du parametre de mou- 
vement au cours du temps dans le champ d 1 observation, on 
calcule au moins un parametre de propagation de 1 • onde de 
cisaillement en au moins certains points du champ d' obser- 
vation ; 
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le parametre de propagation de 1 1 onde de ci- 
saillement qui est calcule au cours de l'etape de carto- 
graphie est choisi parmi la vitesse des ondes de cisaille- 
ment, ie module de cisaillement , 1 ' at tenuat ion des ondes 
5 de cisaillement, l'elasticite de cisaillement, et la vis- 
cosite de cisaillement. 

Par ailleurs, 1* invention a egalement pour objet 
un dispositif d'imagerie utilisant les ondes de cisaille- 
ment pour observer un milieu viscoelast ique diffusant qui 

10 contient des particules ref lechissant les ondes ultrasono- 
res de compression, ce dispositif comprenant des moyens 
d' excitation pour generer une onde elastique de cisaille- 
ment dans le milieu viscoelast ique et des moyens d 1 acqui- 
sition pour observer, au moyen d'au moins une onde ultra - 

15 sonore de compression, le deplacement du milieu visccelas- 
tique soumis a ladite onde de cisaillement, 

caraccerise en ce que les moyens d' excitation sont adaptes 
pour appliquer au milieu viscoelast ique une excitation 
ayant la forme d'une impulsion basse frequence qui pre- 
20 sente une frequence centrale f comprise entre 20 et 5000 
Hz, cecte impulsion basse frequence presentant une duree 
comprise entre l/2f et 20/f, 

en ce que les moyens d 1 acquisition sont adaptes pour ob- 
server la propagation de 1 1 onde de cisaillement simultane- 

25 ment en une multitude de points dans le milieu observe, 
ces points formant un champ d 1 observat ion sensiblement 
continu s'etendant au moins selon un premier axe, lesdits 
moyens d * acquisition etant adaptes pour : 

emettre dans le milieu observe une succession 

30 d'au moins 10 tirs d' ondes ultrasonores de compression a 
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une cadence comprise entre 100 et 100 000 tirs par se- 
conde , 

detecter et enregistrer en temps reel les echcs 
generes par les particules ref lechissantes du milieu vis- 
coelas:ique a chaque tir d'onde ultrasonore, ces echos 
correspondant (directement ou indirectement) a des images 
successives du milieu observe, 

et en ce que ledit dispositif comporte en outre des mcyens 
de traitement d 1 image adaptes pour traiter en temps dif fe- 
re les. . images obtenues a partir des moyens d ' observation , 
au rnoins par intercorrelat ion entre images successives, 
pour determiner en chaque point du champ d ' observat ion un 
parametre de mouvement choisi entre le deplacement et la 
deformation du milieu viscoelast ique , de fagon a obtenir 
ainsi une succession d' images montrant 1' evolution du pa- 
rametre de mouvement du milieu viscoelast ique sous I'effet 
de la propagation de 1 1 onde de cisaillement . 

Dans des modes de realisation preferes du disposi- 
tif selcn 1' invention, on peut eventuellement avoir re - 
cours en outre a 1 1 une et/ou a 1' autre des dispositions 
sui varices : 

les moyens d ' observation comprennent une batte- 
rie de t ransduc teurs qui inclut au moins un transducteur 
et qui est adaptee pour Stre disposee au contact d'une 
surface exterieure delimitant le milieu viscoelast ique , 
les moyens d' excitation etant adaptes pour imposer un de- 
placement impulsionnel £ ladite batterie de transduc- 
teurs ; 

ledit parametre de mouvement est la deformation 
du milieu viscoelast ique ; 
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le champ d 1 observation s'etend au moins selon 
un plan comprenant d'une part, le premier axe et d' autre 
part, un deuxieme axe perpendiculaire au premier axe, la 
batterie de t ransducteurs comportant plusieurs transduc- 
teurs disposes au moins selon le deuxieme axe, des moyens 
de comnande etant prevus pour select ivement faire fonc- 
tionner le dispositif soit dans le mode d'imagerie par on- 
des de cisaillement , soit dans un mode d ' echographie stan- 
dard permettant d'acquerir entre 10 et 100 images par se- 
conde . 

D'autres caracteristiques et avantages de 1' inven- 
tion apparaitront au cours de la description suivante de 
plusieurs de ses formes de realisation, donnees a titre 
d'exemples non limitatifs, en regard des dessins joints. 

Sur les dessins : 

la figure 1 est une vue schematique d'un dispo- 
sitif d'imagerie par ondes de cisaillement selon une forme 
de realisation de 1' invention, 

et la figure 2 est une vue de detail montrant 
une variante du dispositif de la figure 1. 

Sur les differentes figures, les memes references 
designent des elements identiques ou similaires. 

La figure 1 represente un example de dispositif 
d'imagerie par ondes de cisaillement selon 1 ' invention, 
pour ecudier la propagation des ondes elastiques de ci- 
saillement dans un milieu viscoelast ique 1 qui est diffu- 
sant vis a vis des ondes ultrasonores de compression, et 
qui peut etre par exemple : 
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un corps inerte, notamment dans le cas du con- 
trole de qualite pour des applications industriel les , no- 
^ammeac agro-al iment aires , 

ou un corps vivant, par exemple une partie du 
5 corps d'un patient, dans le cas des applications medica- 
les . 

Ce dispositif comporte un transducteur acoustique 
tel qu'un haut-parleur 2 ou un pot vibreur qui est dispcse 
contra la surface extSrieure 3 du milieu observe 1, cette 
10 surface 3 etant constitute par example par la peau du'pa- 
tient dans les applications medicales.. 

Le haut-parleur 2 peut etre commande par un micro- 
ordinateur 4, par exemple par 1 1 intermediaire d'un circuit 
generateur d r impulsions basse frequence G (ce circuit peut 

15 etre const itue notamment par la carte son du micrc- 
ordinateur 4} et d'un ampl if icateur A, pour appliquer a la 
surface 3 du milieu observe une excitation sous la forme 
d'une impulsion basse frequence, de facpon a generer une 
onde de cisaillement dans le milieu viscoelast ique 1. 

-0 Cette impulsion basse frequence presente en gene- 

ral une amplitude pouvant etre de 1 ' ordre de 1 mm et une 
frequence centrale f comprise entre 20 et 5000 Hz, appli- 
quee pendant une duree comprise entre l/2f et 20/f. De 
preference, la duree d 1 application de 1 ' impulsion basse 

25 frequence est comprise entre l/2f et 2f et la frequence f 
est comprise entre 30 et 1000 Hz, cette frequence pcuvant 
typiquement etre de 1' ordre de 50 Hz. 

En variante 1 1 onde acoustique de cisaillement 
pourrait egalement etre obtenue (avec les caracterist iques 

30 susmentionnees d» amplitude et de frequence) 
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par un choc genere localement par contact phy- 
sique au niveau de la surface exterieure du milieu visco- 
elastique, par au moins un actionneur mecanique automati- 
que commande par le micro-ordinateur 4, 

par un choc applique manuellement par contact 
physique au niveau de la surface exterieure du milieu vis- 
coelast ique , 

ou encore, par un mouvement impulsionnel natu- 
ral d'un organe interne du corps humain ou animal (par 
example un battement du coeur) dans les applications medi- 
cales . 

L'onde elastique de cisaillement produite par le 
haut-parleur 2 se deplace avec une vitesse Cs relativement 
faibie, de 1 1 ordre de quelques m/s ( typiquement , de 1 a 10 
m/s dans le corps humain) en produisant des mouvements in- 
ternes dans le milieu viscoelastique 1 observe. 

Ces mouvements sont suivis en envoyant dans le mi- 
lieu 1, des ondes ultrasonores de compression qui inter- 
agissenc avec les particules diffusantes 5 contenues dans 
le milieu 1, lesquelles particules sont ref lechissantes 
pour les ondes ultrasonores de compression. Les particules 
5 peuvent etre constitutes par toute heterogenei te du mi- 
lieu 1, et notamment, lorsqu'il s 1 agit d'une application 
medicale, par des particules de collagene presentes dans 
les tissus humains. 

Pour observer la propagation de l'onde de ci- 
saillement, on utilise done une sonde ultrasonore 6 dispo- 
see contre la surface exterieure 3 du milieu observe 1. 
Cette sonde envoie, selon un axe X, des impulsions d* ondes 
ultrasonores de compression du type de celles couramment 
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utilisees en echographie, a une frequence comprise par 
exemple entre 1 et 100 MHz et de preference entre 3 et 15 
MHz. On notera que la sonde 6 peut etre disposee : 

soit du meme c6te du milieu 1 que le haut- 
parleur 2, comme represents sur la figure 1, 

soit at 1' oppose du haut-parleur 2 par rapport 
au milieu 1, 

soit dans toute autre position, par exemple 
dans une disposition transversale par rapport au haut- 
parleur 2. — 

La sonde ultrasonore 6 est constitute par une bat- 
terie de n transducteurs ultrasonores Tl, T2 , Ti, 
. .., Tn, n etant un nombre entier n au moins egal a 1. 

Cette sonde 6 se presente le plus souvent sous la 
forme d'une barret te linSaire pouvant comprendre par exem- 
pt n = 128 transducteurs alignes selon un axe Y perpendi- 
culaire a I'axe X, envoyant simul tanement leurs impulsions 
d'ondes ultrasonores de fa<;on a generer une onde "plane" 
(c'est a dire en 1* occurrence une onde dont le front 
d 1 onde est rectiligne dans le plan X, Y) ou tout autre 
type d'onde eclairant 1' ensemble du champ d ' observat ion . 

En variante, la batterie 2 de transducteurs peut 
eventuellement etre reduite a un transducteur unique Tl, 
ou au contraire se presenter sous la forme d 1 un reseau a 
deux dimensions s'etendant par exemple selon un plan per- 
pendiculaire a 1 1 axe X. 

Chacun des transducteurs Tl, T2 , ... Tn est corn- 
mande par le micro-ordinateur 4 ou par une unite centrale 
CPU {qui est contenue par exemple dans une bale electronic 
que 7 reliee par un cable souple a la sonde 6), pour erne:- 
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tre des tirs successifs d'ondes de compression ultrascno- 
res dans le milieu 2, au cours d'une phase d'observation 
qui peut dur.er par exemple moins d'une seconde et au cours 
de laquelle sont emis p tirs d'ondes ultrasonores de com- 
pression (p etant un entier compris entre 100 et 10 000 et 
de preference entre 1000 et 100 000) , a une cadence com- 
prise entre 100 et 100 000 tirs par seconde et de prefe- 
rence comprise entre 1000 et 100 000 tirs par seconde, no- 
tamment entre 1000 et 10 000 tirs par seconde (cette ca- 
dence est limitee par le temps d' aller-retour de I'onde de 
compression dans le milieu 1, done par l'epaisseur du mi- 
lieu 1 dans la direction X : il faut en effet que tous les 
echos generes par 1 » onde de compression aient ete rectus 
par la sonde 6 avant d 1 envoy er une nouvelle onde de com- 
pression) . 

Les tirs d'ondes ultrasonores de compression de la 
phase d'observation commencent de preference juste avant 
1' emission de 1 ' onde de cisaillement . 

De plus, dans le cas ou 1 ' onde de cisaillement est 
generee par un mouvement impulsionnel d'un organe d'un 
corps vivant, on peut avantageusement synchroniser le de- 
marrage des tirs d'ondes ultrasonores de compression avec 
ce mouvement impulsionnel. par exemple dans le cas ou 
1'onde de cisaillement est generee par un battement de 
caeur, on peut synchroniser le demarrage des tirs d'ondes 
ultrasonores de compression avec une phase choisie de 
1 ' electrocardiogramme . 

Chacun de ces tirs donne lieu a la propagation 
d'une onde ultrasonore de compression dans le milieu 1, 
avec une vitesse de propagation beaucoup plus elevee que 
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les ondes de cisaillement , par exemple de 1 1 ordre de 1500 
m/s dans le corps humain. 

L'onde ultrasonore ainsi generee interagit avec 
les particules ref lechissantes 5, ce qui genere des echos 
ou aucres perturbations analogues du signal, connus en soi 
sous le nom de "bruits de speckle" dans le domaine de 
1 * echographie . 

Ces "bruits de speckle" sont captes par les trans- 
ducteurs Tl, ... ,Tn aprds chaque tir. Le signal sij (t) ain- 
si capte par chaque transducteur Ti apr§s. le tir n° j est 
tout d'abord echant illonne £ haute frequence (par exemple 
de 30 a 100 MHz) et numerise en temps reel (par exemple 
sur 8 bits, ou dans certains cas sur 1 bit) par un echan- 
tillonneur appartenant a la baie 7 et reli£ a ce transduc- 
teur, respect ivement El, E2, ... En. 

Le signal sij(t) ainsi echant illonne et numerise 
est ensuite memorise, egalement en temps reel, dans une 
memoire Mi appartenant 3 la baie 7 et propre au transduc- 
teur Ti . 

Chaque memoire Mi presente par exemple une capaci- 
te de I 1 ordre de 1 Mo, et contient 1' ensemble des signaux 
sij (t) regus successivement pour les tirs j = 1 a p. 

En temps differe, apres la memorisation de tous 
les signaux sij (t) correspondant a une meme propagation 
d'onde de cisaillement, 1' unite centrale CPU fait retrai- 
ter ces signaux par un circuit sommateur S appartenant a 
la baie 7 (ou bien elle effectue elle-meme ce traitement, 
ou encore ledit traitement peut etre effectue dans le mi- 
cro-ordinateur 4) , par un processus classique de formation 
de voies. 
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On genere ainsi des signaux Sj (x, y) qui corres- 
pondent chacun a 1' image du champ d' observation apres le 
tir n° j . 

Par exemple, on peut determiner un signal Sj (t) 
5 par la formule suivante : 

SAO = 2>. ( x > yy si JV(x, y) + d t (x, v) / v] 

OU : 

sij est le signal brut pergu par le transduc- 
teur n° i apres le tir d'onde ultrasonore de compression 
10 n° j , 

t(x,y) est le temps mis par 1 ' onde ultrasonore 
de compression pour atteindre le point du champ 
d' observation de coordonnees (x,y), avec t = 0 au debut du 
tir n° j , 

15 - di(x,y) est la distance entre le point du champ 

d' observation de coordonnees (x,y) et le transducteur n° 
i, ou une approximation de cette distance, 

V est la vitesse moyenne de propagation des on- 
des acoustiques ultrasonores de compression dans le milieu 
20 viscoelastique observe. 

et ai(x,y) est un coefficient de ponderation 
tenant compte de lois d' apodisation (en pratique, on pour- 
ra dans de nombreux cas considerer que ai(x,y)=l) . 

La formule ci-dessus s' applique mutatis mutandis 
25 lorsque le champ d' observation est a 3 dimensions (reseau 
plan de transducteurs) , en remplac^ant les coordonnees spa- 
• tiales (x,y) par (x,y,z). 
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Lorsque la sonde 6 comport e un seul transducteur , 
l'etape de formation de voies est inutile, et on a direc- 
tement Sj (x) = sj[2.x/V], avec les memes notations que 
ci -dessus . 

5 Apres l'etape eventuelle de formation de voies, 

1' unite centrale CPU memorise dans une memoire centrale M 
apparcenant a la baie 7 les signaux d' images Sj(x,y) ou 
Sj (x) ou Sj(x,y,z), qui correspondent chacun au tir n° j. 
Ces signaux peuvent egalement etre memorises dans le mi- 
10 cro-ordinateur 4 lorsqu'il effectue lui-meme le traitement 
d' image . 

Ces images sont ensuite traitees deux a deux, ou- 
jours en temps diff§re, par intercorrelat ion . Cette inter- 
correlation peut etre realisee dans un circuit DSP appar- 

15 tenant a la baie 7, ou etre programmee dans 1' unite cen- 
trale CPU ou dans le micro-ordinateur 4. 

A titre d ! exemple, cette intercorrelat ion peut se 
faire en comparant les signaux Sj (x,y) et Sj+1 (x,y) 
(dans le cas d'un champ d' observation a 2 dimensions) sur 

20 des fenetres spatiales glissantes de longueur Ax predeter- 
minee, pouvant aller par exemple de X a 10X. , oQ X est la 
longueur d'onde des ondes ultrasonores de compression 
(soit environ 0,42 S 4,2 mm a 3,5 MHz dans l'eau ou le 
corps humain) . De plus, les fenetres susment ionnees peu- 

25 vent se recouvrir l'une 1 1 autre sur environ 20% de leur 
longueur selon 1 ' axe X. 

Au cours de ce processus d' intercorrelat ion, on 
maximise une fonction d ' intercorrelat ion 

<Sj (x,y) ,Sj+l (x,y) > afin de determiner le deplacement subi 
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par chaque particule 5 dormant lieu a un echo ultrasonore , 
dans la direction X. 

Des exemples de tels calculs d' intercorrelat ion 
sont donnes dans 1'etat de la technique, notamment par 
5 O'Donnell et al . (« Internal displacement and strain ima- 
ging using speckle tracking », IEEE transactions on ultra- 
sonic, f erroelectrics, and frequency control, vol. 41, n° 
3, mai 1994, p. 314-325) et par Ophir et al . 
(« Elastography : a quantitative method for imaging the 

10 elasticity of biological tissues », Ultrasonic imag., vol. 
13 , p. 111-134, 1991) . 

On obtient ainsi une succession de champs de de- 
placements Djx(x,y) du milieu 1 sous l'effet de 1 ' onde de 
cisaillement , dans la direction X. 

15 Cette succession de champs de displacements est 

stockee dans la memoire M ou dans le micro-ordinateur 4 et 
peut etre visualisee, notamment au moyen de l'ecran 4a du 
micro-ordinateur, sous la forme d'un film ralenti ou la 
valeur des deplacements est illustree par un parametre op- 

20 tique tel que par un niveau de gris ou par un niveau chro- 
mat ique . 

On visualise ainsi parfaitement les differences de 
propagation de 1 • onde de cisaillement entre les zones de 
caracteristiques differentes du milieu 1, par exemple les 
25 tissus sains et les tissus cancereux dans le cas d'une ap- 
plication medicale. 

Ce film de propagation de 1 ■ onde de cisaillement 
est en outre superposable avec une image echographique 
classique, qui peut etre generee par le dispositif decrit 
30 ci-dessus, lequel est apte a fonctionner : 
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soit dans le mode d'imagerie par ondes de ci- 

saillement , 

soit dans un mode d ' echographie standard, en 
fonction des commandes re?ues par exemple depuis le cla- 
5 vier 4b du micro-ordinateur . 

Par ailleurs, il est egalement possible de calcu- 
ler non pas les d£placements de chaque point du milieu ob- 
serve 1, mais les deformations Ejx(x,y) du milieu dans la 
direction X, c'est a dire les derivees des deplacements 

10 Djx(x,y) par rapport & x. 

Ces champs de deformations successifs sont utili- 
sables comme precedemment pour visualiser clairement la 
propagation de 1 ' onde de cisaillement sous la forme d 1 un 
film, et presentent en outre 1'avantage de s'affranchir 
15 des deplacements de la sonde 6 par rapport au milieu ob- 
serve 1 . 

Cette variante est part iculierement interessante 
dans la forme de realisation de la figure 2, ou la sonde 6 
est portee par le haut-parleur ou vibreur 2, ce qui impli- 

20 que forcement des mouvements de ladite sonde puisque c'est 
elle-meme qui gendre alors 1 1 onde de cisaillement. 

A partir des champs de deplacements ou de deforma- 
tions, on peut le cas echeant proceder en outre a une 
etape de cartographie au cours de laquelle, a partir de 

25 1' evolution du parametre de mouvement {deplacement ou de- 
formation) au cours du temps dans le champ d 1 observation 
X, Y (ou X dans la ces d'un seul transducteur , ou X, Y, Z 
dans le cas d'un reseau plan de transducteurs) , on calcule 
au moins un parametre de propagation de l'onde de ci- 

30 saillement, soit en certains points du champ d ' observation 
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choisis par 1 1 ut il isateur a partir du micro-ordinateur 4, 
soit dans tout le champ d 1 observation . 

Le parametre de propagation de l'onde de cisaille- 
ment qui est calcule au cours de l'etape de cartographie 
5 est choisi par exemple parmi la vitesse Cs des ondes de 
cisaillement, le module de cisaillement yi, 1 ' attenuation a 
des ondes de cisaillement, l'elasticite de cisaillement 
ILilviscosite de cisaillement \x2 . 

Ce calcul est effectue par un processus classique 
10 d* inversion, dont un exemple est donne ci-apres dans le 
cas d'un champ d 1 observation a deux dimensions (le meme 
processus s 1 appliquerait mutatis mutandis dans le cas d'un 
champ d' observation a une ou a trois dimensions, respec- 
tivement pour un seul transducteur Tl ou pour un reseau 
15 plan de transducteurs) . 

Dans cet exemple on fera 1 1 approximation que la 
viscosite de cisaillement \x2 est nulle et que le milieu 
est isotrope. 

L' equation d'onde donnant le vecteur deplacement D 
20 de chaque point du milieu 1 s'ecrit : 

d z D 

/7— = A(//.D) (I), 

Ct ' 

ou p est la densite du milieu 1, yi est le module de ci- 
saillement (par hypothese reduit a sa partie reelle 
l'elasticite de cisaillement ^1 puisque la viscosite de 
25 cisaillement (i2 est supposee nulle) . 

Pour la premiere composante u du vecteur D, c'est 
a dire pour le deplacement du milieu 1 dans la direction 
X, on a done : 
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d 2 u d 2 (/ju) d 2 (yu) 
df dx~ dy 2 

Apres transformee de Fourier temporelle et discre- 
tisation, cette equation peut etre ecrite sous la forme 
matricielle suivante, qui peut etre ecrite pour chaque 
5 frequence du spectre de l'onde de cisaillement : 

B = HM (III) , ou : 

M est un vecteur de dimension (L+2 ) . (M+2 ) -4 
dont chaque composante vaut |i lm , c'est a dire la valeur 
locale~"du module de cisaillement en chaque point discreti- 

10 se de coordonn§es (xi,y m ) , oil 1 et m sont des entiers com- 
pris entre 0 et respect ivement L+l et M+l en eliminant les 
couples (l,m) valant (0,0), (0,M+1), (L+1,0) et (L+l, M+l), 
L+2 et M+2 etant les nombres de points discretises dans 
1' image du milieu 1 respect ivement selon les axes X et Y, 

15 B est un vecteur de dimension L.M, dont les 

composantes valent* -tn 2 .p.Ui m ou w est la pulsation de 
l'onde de cisaillement basse frequence, p est la densite 
du milieu, Ui^ est la transformee de Fourier temporelle du 
deplacement u au point de coordonnees (xi,y m ) , 1 etant 

20 compris entre 1 et L et m etant compris entre 1 et M, 

et H est une matrice de dimension L.M lignes 
sur (L+2) . (M+2) -4 colonnes, dont les composantes sont ega- 
lement toutes connues £ partir de 1' equation d'onde. 

En juxtaposant un nombre suffisant d' equations 

25 (III) correspondant respectivement a differentes disponi- 
bles du spectre de frequences de l'onde de cisaillement, 
on obtient une equation matricielle globale qui peut etre 
resolue par inversion matricielle, pour obtenir le vecteur 
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M, c'est a dire la valeur du module de cisaillement ji en 
tout point du champ d' observation . 

On peut alors en deduire eventuellement la valeur 
locale de la vitesse de propagation Cs de 1 ' onde de ci- 
saillement en chaque point, par la formule : 




Le mode de calcul serait le meme en utilisant non 
plus les deplacements , mais les deformations du milieu ob- 
serve 1 . 
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RFVTTMDTCATIONS 

1. Procede d'imagerie utilisant les ondes de ci- 
saillement pour observer un milieu viscoelastique diffu- 
5 sant (1) qui contient des particules (5) ref lechissant les 
ondes ultrasonores de compression, procede dans lequel on 
genere une onde elastique de cisaillement dans le milieu 
viscoelastique et on observe, au moyen d'au moins une onde 
ultrasonore de compression, le deplacement du milieu vis- 
it) coelastique (1) soumis & ladite onde de cisaillement, 

caractirise en ce qu'on genere 1 1 onde de cisaillement en 
appliquant au milieu viscoelastique une excitation ayant 
la forme d'une impulsion basse frequence qui presente une 
frequence centrale f comprise entre 20 et 5000 Hz, cette 
15 impulsion basse frequence presentant une duree comprise 
entre l/2f et 20/f, 

en ce qu'il comporte une etape d ' observation de propaga- 
tion au cours de laquelle on observe la propagation de 
l'onde de cisaillement simultanement en une multitude de 

20 points dans le milieu observe, ces points formant un champ 
d 1 observation sensiblement continu s'etendant au moins se- 
lon un premier axe (X), cette etape d • observation de pro- 
pagation de l'onde de cisaillement consistant a : 

emettre dans le milieu observe une succession 

25 d'au moins 10 tirs d 1 ondes ultrasonores de compression a 
une cadence comprise entre 100 et 100 000 tirs par se- 
conde, 

detecter et enregistrer en temps reel les echos 
generes par les particules ref lechissantes du milieu vis- 
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coelastique a chaque tir d'onde ultrasonore, ces echos 
correspondant a des images successives du milieu observe, 
et en ce que ledit precede comporte en outre une etape ul- 
terieure de traitement d' image au cours de laquelle on 
traite en temps differe les images ainsi obtenues au moins 
par intercorrelation entre images successives, pour deter- 
miner en chaque point du champ d' observation un parametre 
de mouvement choisi entre le deplacement et la deformation 
du milieu viscoelast ique , de fagon a obtenir ainsi une 
succession d' images montrant 1 'evolution du parametre de 
mouvement du milieu viscoelast ique sous l'effet de la pro- 
pagation de I'onde de cisaillement . 

2. Procede selon la revendication 1, dans lequel 
la duree de 1' impulsion basse frequence est comprise entre 
l/2f et 2f. 

3. Procede selon la revendication 1 ou la revendi- 
cation 2, dans lequel la frequence centrale de 1' impulsion 
basse frequence est comprise entre 30 et 1000 Hz. 

4. Procede selon 1 ' une quelconque des revendica- 
tions precedentes, dans lequel le milieu viscoelast ique 
(1) observe est const itue par un corps vivant comprenant 
au moins un organe interne soumis a des mouvements impul- 
sionnels, 1' impulsion basse frequence qui genere I'onde de 
cisaillement etant constitute par un mouvement impulsion- 
nel dudit organe interne. 

5. Procede selon 1 1 une quelconque des revendica- 
tions 1 a 3, dans lequel le milieu viscoelast ique (1) ob- 
serve est delimite par une surface exterieure (3) , et 

*1' impulsion basse frequence est appliquee au niveau de 
cette surface exterieure. 
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6. Proc£de selon la revendication 5, dans lequel 
1' impulsion basse frequence f*st appliquee par un mcyen 
d' excitation choisi parmi : 

une onde acoustique generee par au moins un 
transducteur acoustique (2) , 

et un choc genere localement par contact physi- 
que au niveau de la surface exterieure (3) du milieu vis- 
coelastique . 

7. Procede selon la revendication 5 ou la revendi- 
cation 6, dans lequel on §met les tirs d'ondes ultrasono- 
res de compression et on detecte les echos generes par les 
particules ref lechissantes du milieu viscoelast ique au 
moyen d 1 une bacterie de transducteurs (6) qui comprend au 
moins un transducteur (Tl, Tn) et qui est disposee au 
contact de la surface exterieure du milieu viscoelast ique , 
1 1 onde de cisaillement etant appliquee au milieu visco- 
elast ique en imposant un deplacement impulsionnel a iadite 
batterie de transducteurs. 

8. Procede selon la revendication 7, dans lequel 
ledit parametre de mouvement est la deformation du milieu 
viscoelastique . 

9. Procede selon 1 1 une quelconque des revendica- 
tions precedentes, dans lequel, au cours de 1 ' observation 
de la propagation de l'onde de ci saillement, on emet entre 
100 et 10 000 tirs d'ondes ultrasonores de compression a 
une cadence comprise entre 1000 et 100 000 tirs par se- 
conde . 

10. Procede selon 1 • une quelconque des revendica- 
tions precedentes, dans lequel le champ d 1 observation 
s'etend au moins selon un plan comprenant d'une part, le 
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premier axe (X) et d' autre part, un deuxieme axe perpendi- 
culaire au premier axe (Y) . 

11. Proced<§ selon la revendication 10, dans le- 
quel, au cours de 1' etape d 1 observation de propagation, on 

5 utilise une batterie (6) de plusieurs transducteurs acous- 
tiques (Tl, Tn) disposes au moins selon le deuxieme axe 
(Y) pour emettre les tirs d'ondes acoustiques ultrasonores 
de compression et detecter les echos generes par les par- 
ticules r<§f lechissantes (5) du milieu viscoelastique, les 

10 echos detectes par chaque transducteur acoustique etant 
memorises directement sans traitement pr<§alable au cours 
de 1 1 etape d ' observation de propagation, et 1 ' etape de 
traitement d f image comprenant une sous-etape preliminaire 
de formation de voies au cours de laquelle on genere une 

15 image du milieu viscoelastique correspondant a chaque tir 
d'onde ultrasonore de compression, par combinaison d T au 
moins certains des echos rectus par les differents trans- 
ducteurs (Tl, Tn) . 

12. Procede selon 1 ' une quelconque des revendica- 
20 tions precedentes, dans lequel 1 1 etape de traitement 

d f images est suivie par une etape de visualisation au 
cours de laquelle on visualise au ralenti un film consti- 
tue par la succession des images traitees, chaque point de 
chaque image presentant un parametre optique qui varie 
25 suivant la valeur du parametre de mouvement affecte a ce 
point . 

13. Procede selon la revendication 12, dans lequel 
le parametre optique est choisi parmi le niveau de gris et 
le niveau chromatique. 
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14. Procede selon 1 1 une quelconque des revendica- 
tions precedentes, dans lequel 1 ' etape de traitement 
d' images est suivie par une etape de cartographie au cours 
de laquelle, a partir de devolution du parametre de mou- 

5 vement au cours du temps dans le champ d ' observation, on 
calcule au moins un parametre de propagation de l'onde de 
cisaillement en au moins certains points du champ d'obser- 
• vat ion . 

15. Proc§de selon la revendication 14, dans lequel 
10 le parametre de propagation de l'onde de cisaillement qui 

est calcule au cours de 1 ' etape de cartographie est choisi 
parmi la vitesse des ondes de cisaillement, le module de 
cisaillement, 1 ' attenuation, des ondes de cisaillement, 
l'elasticite de cisaillement, et la viscosite de cisaille- 
15 ment . 

16. Dispositif d'imagerie utilisant les ondes de 
cisaillement pour observer un milieu viscoelast ique diffu- 
sant (1) qui contient des particules (5) ref lechissant les 
ondes ultrasonores de compression, ce dispositif compre- 

20 nant des moyens d f excitation (2) pour generer une onde 
elastique de cisaillement dans le milieu viscoelast ique et 
des moyens d 1 acquisition (CPU, Ti, Ei, Mi) pour observer, 
au moyen d 1 au moins une onde ultrasonore de compression, 
le deplacement du milieu viscoelast ique soumis a ladite 

25 onde de cisaillement, 

caracterise en ce que les moyens d' excitation (2) sont 
adaptes pour appliquer au milieu viscoelast ique (1) une 
excitation ayant la forme d'une impulsion basse frequence 
qui presente une frequence centrale f comprise entre 20 et 
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5000 Hz, cette impulsion basse frequence presentant une 
duree comprise entre l/2f et 20/f, 

en ce que les moyens d 1 acquisition (CPU, Ti, Ei, Mi) sont 
adaptes pour observer la propagation de 1 ' onde de ci- 
5 saillement simultanement en une multitude de points dans 
le milieu observe, ces points formant un champ d 1 observa- 
tion sensiblement continu s'etendant au moins selon un 
premier axe (X), lesdits moyens d r acquisition etant adap- 
tes pour : 

10 - emettre dans le milieu observe une succession 

d'au moins 10 tirs d'ondes ultrasonores de compression si 
une cadence comprise entre 100 et 100 000 tirs par se- 
conde , 

detecter et enregistrer en temps r£el les echos 
15 generes par les particules ref lechissantes du milieu vis- 
coelastique 3. chaque tir d'onde ultrasonore, ces #chos 
correspondant a des images successives du milieu observe, 
et en ce que ledit dispositif comporte en outre des moyens 
de traitement d f image (CPU, S, DSP) adaptes pour traiter 
20 en temps differe les images obtenues a partir des moyens 
d ' observation , au moins par intercorrelat ion entre images 
successives, pour determiner en chaque point du champ 
d' observation un paramdtre de mouvement choisi entre le 
deplacement et la deformation du milieu viscoelast ique , de 
25 fa^on a obtenir ainsi une succession d' images montrant 
1' evolution du parametre de mouvement du milieu viscoelas- 
tique sous I'effet de la propagation de l f onde de ci- 
saillement . 

17. Dispositif selon la revendicat ion 16, dans le- 
30 quel les moyens d 1 observation comprennent une batterie (6) 
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de transducteurs (Tl, Tn) qui inclut au moins un trans- 
ducteur et qui est adaptee pour ecre disposee au contact 
d'une surface ext£rieure (3) d<§limitant le milieu visco- 
elastique (1), les moyens d'excitation (2) etant adaptes 
pour imposer un deplacement impulsionnel a ladite batterie 
(6) de transducteurs. 

18. Dispositif selon la revendicat ion 17, dans le- 
quel ledit parametre de mouvement est la deformation du 
milieu visco£lastique . 

19. Dispositif selon 1 1 une quelconque des revendi- 
cations 17 et 18, dans lequel le champ d ■ observation 
s'etend au moins selon un plan comprenant d'une part, le 
premier axe (X) et d* autre part, un deuxieme axe (Y) per- 
pendiculaire au premier axe, la batterie (6) de transduc- 
teurs comportant plusieurs transducteurs (Tl, Tn) dis- 
poses au moins selon le deuxieme axe, des moyens de com- 
mande (4, CPU) etant pr6vus pour select ivement faire fonc- 
tionner le dispositif soit dans le mode d'imagerie par on- 
des de cisaillement , soit dans un mode d ' echographie stan- 
dard permettant d'acqu^rir entre 10 et 100 images par se- 
conde . 
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